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When a very thin capiUary is inserted into a liquid, the liquid is sucked into it: this imbibition 
process is controlled by a balance of capillary and drag forces, which are hard to quantify experimen- 
tally, in particularly considering flow on the nanoscale. By computer experiments using a generic 
coarse-grained model, it is shown that an analysis of imbibition forced by a controllable external 
pressure quantifies relevant physical parameter such as the Laplace pressure, Darcy's permeability, 
effective pore radius, effective viscosity, dynamic contact angle and slip length of the fluid flowing 
into the pore. In determining all these parameters independently, the consistency of our analysis of 
such forced imbibition processes is demonstrated. 



Flowing fluids confined to pores with diametera on the 
fim to nm scale are important in many contexts: oil 
recovery from porous rocks separation processes in 
zeolithes Q; nanofluidic devices such as hquids in nan- 
otubes [1]; nanohthography [3] jUanolubrication fluid 
transport in living organisms [6| and many other appli- 
cations [H- However, despite its importance for so many 
processes in physics, chemistry, biology and technology, 
the flow of fluids into (and inside) nanoporous materials 
often is not well understood: the effect of pore surface 
structure on the flowing fluid is difficult to assess, 

in terms of hydrodynamics, the problem differs dramat- 
ically from the macroscopic fiuid dynamics 0, [l^ ; and 
although very beautiful experiments have recently been 
made (e.g. [ll El), 

more information is needed for a 
complete description of the relevant microfiuidic process. 

In the present work, we propose to use forced imbibition 
with the external pressure as a convenient control param- 
eter to obtain a much more diverse information on the 
parameters controlling flow into capillaries than hereto- 
fore possible. Extending our recent study of imbibition 
at zero pressure , we concisely describe the theoretical 
basis for this new concept, and provide a comprehensive 
test of the concept in terms of a " computer experiment" 
on a generic model system (a fluid composed of Lennard- 
Jones particles flowing into a tube with a perfectly crys- 
talline (almost rigid [iSj]) wall and spherical cross-section, 
see Fig. [1] We also provide a stringent test of our descrip- 
tion by having estimated all parameters of the theory in 
independent earlier work [l3|), and hence there are no 
adjustable parameters whatsoever. We emphasize that 
our procedures and analysis could be followed in experi- 
ments with real materials fully analogously. 

We briefly summarize the pertinent theory. On a 
macroscopic scale the rise of a fluid meniscus at height 
H(t) over the entrance of a capillary with time t is de- 
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FIG. 1: A snapshot of the capillary during imbibition (cross 
section along the axis of the cylindrical capillary oriented in 
a;-direction) . At the left side there is a box, which initially has 
a cubic shape with linear dimensions of 40 Lennard- Jones di- 
ameters a, and contains A'' = 40000 particles interacting with 
the Lennard- Jones potential ULj{r) = 4e[((T/r)^'^ — (a/r)^], 
r being the distance between the particles, e — lA (in tem- 
perature units, ksT = 1). The left wall of the box can move 
to maintain a given pressure P^xt in the reservoir box. The 
capillary has a radius R — 11a and length SOcr. It ends at 
a hard wall on the right side while it is open toward the left 
side, to connect with the reservoir box, which has a circular 
opening in the right wall (shown in yellow). 



scribed (at zero external pressure) by the well-known 
Lucas- Washburn equation 



cos 9 



(1) 



Here jlv is the liquid-vapour surface tension of the liq- 
uid, T] the shear viscosity of the fluid, R the pore ra- 
dius, 9 the contact angle, and Hq a constant (which ac- 
counts for the fact that Eq.([T]) holds only after some 
transient time when inertial effects have already van- 
ished). Eq.(IT]) follows when one balances the viscous 
drag force ^H(t) '^^Jf'' with the Laplace pressure Pl = 
2-fLV cos{9)/R. 
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The applicability of Ec^ffl for ultrathin pores has been 
rather controversial 

0,1 

ItI. [l8|. This debate was clari- 
fied [ l 3| by recalling that on the nanoscale the slip length 
S [l^, 120] must not be neglected. According to the defi- 
nition of this length, the drag force under slip flow con- 
ditions in a tube of radius R and slip length S is equal 
to the drag force for a no-slip flow in a tube of effective 
radius R + S. Thus one ends up with a modified Lucas- 
Washburn relationship: 
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On the left hand side of Eq. [5] we have now also included 
an external pressure term Pext- If one uses Darcy's per- 
meability )^ K = [R + 5f/S, Eq. [2] can be written in 
a form which does not depend on the capillary radius 

anymore, introducing also the rate v{Pext) — ~~dt~^j 
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For constant Pext, Eq. [2]is easily integrated to 
2k , 



H\t) = —{PL+Pext)t + Hl 



(3) 
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Eq. [3] shows that v{Pext) varies linearly with Pext, so 
measuring this relationship yields both parameters Pl 
and r]/K. Instead of using the height H{t) by observ- 
ing the meniscus, one may alternatively estimate v{Pext) 
from the time variation of the mass of the fluid inside the 
capillary (i.e. the total number of particles N{t) oc H{t) 
which has entered the capillary up to the time t). In 
contrast, the classical experiments on spontaneous imbi- 
bition of a fluid, where Pext = 0, yield only the product 
kPl, and hence even if the fluid viscosity r] is known, 
one cannot discern the effects due to the driving force 
(oc Pl) and due to the drag force (oc k). Moreover, Eq.[3] 
suggests the intriguing possibility of applying the present 
concepts to the most general case of porous media [ll , ir- 
respective of the particular geometry and topology of the 
channels in such materials, but this will not be followed 
up here. 

We now present a test of the above concepts by a 
quantitative analysis of the computer experiment out- 
lined in Fig. [T] We assume also a Lennard- Jones inter- 
action (of strength ewL between the wall and the fluid 
particles (see (l3[ for details on how the wall is atomisti- 
cally modeled), and study the cases of both nonwettable 
i^WL < 0.65) and wettable {ewL > 0.65) walls. Fig. [2] 
shows a plot of v{Pext) vs. Pext- One can see that there 
is a broad regime where the variation of v{Pext) with 
Pext is indeed linear (deviations from linearity for large 
Pext can be attributed to a slight increase of viscosity 
with increasing fluid density at large pressures). Thus 
Fig. [2] demonstrates that indeed a rather precise estima- 
tion of both Pl and k is possible. This is important in 
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FIG. 2; Imbibition rate v — dH^{t)/dt plotted vs. exter- 
nal pressure Pe.xt, for two choices of ewL (as indicated) that 
controls the wettability of the capillary wall: for ewL = 0.4 
the contact angle cos 9 ~ —0.56 (the fluid does not wet the 
wall); for ewL = 1.2 one finds cos9 ~ 0.92, i.e. almost com- 
plete wetting [l^. The dashed lines show fits to the observed 
v{Pe.xt) to straight lines, reducing the range to Pext < 0.6. 
The slope of the straight lines is equal to ^k/t], the intersec- 
tion point with the axis v = yields the Laplace pressure 
Pl- Note that the Molecular Dynamics simulation applies 
the velocity- Verlet algorithm with a timestep 5t = 0.01/\/48, 
choosing the particle mass m = 1 and units where a = 1 
and fesT = 1. 



many cases, e.g. nanocapillaries or porous media, where 
neither Pl nor k can be reliably predicted theoretically 
(because information is missing, e.g. the effective chan- 
nel radius R or the (dynamic) contact angle 6 or the slip 
length 5 may be unknown). 

From the Laplace pressure Pl one can readily obtain 
information on the contact angle (if interfacial tension 
^LV and pore radius R are known). Fig. [3] shows the 
variation of Pl with ewL in our model. By "measur- 
ing" the contact angle 9 dependence on twL in a sep- 
arate simulation, as well as 7lv, we can predict Pl as 
Pl = 2jLvcos6/ R, as noted above. Fig.[3]shows that the 
agreement between this prediction and the observations 
is excellent. 

Fig. [4] shows that also the "friction coefficient" of the 
imbibition (per unit length) ri/2K strongly depends on 
the wettability of the pore wall. The computer exper- 
iment has the bonus that it yields insight into the be- 
havior of the system on the nanoscale in arbitrary de- 
tail. This is demonstrated by the density profiles of the 
moving meniscus, shown for ewL = 0.4 and ewL = 1-2, 
respectively. While no layering of the fluid is observed 
in the case of nonwetting fluids, ewL = 0.4, for a wet- 
table wall the profile for ewL = 1.2 indicates significant 
density oscillations in the vicinity of the wall, i.e. fluid 
"layering" 

Since the shear viscosity 7? has been determined inde- 
pendently for our system 13|, rj = 6.34 ± 0.15, the ratio 
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FIG. 3: Plot of the Laplace pressure Ph against the strength 
^WL of the wall-liquid interaction. Triangles denote estimates 
from the meniscus position H{t), squares are derived from the 
number of particles N{t) that have entered the capillary. Full 
circles and a dashed line show the theoretical prediction of 
Pl = 2-ycosd/R, cf. text. The estimate 'jlv = 0.735 ± 0.015 
was taken from [13l |. 



r]/2K is readily converted into an estimate for the slip 
length 6 (Fig.HjD). The gradual decrease of S with grow- 
ing wettability of the wall ewL is clearly demonstrated. 

In conclusion, we have modelled a possible and simple 
experimental set-up (Fig.[T]) by computer simulation and 
provided a theoretical framework, by slightly extending 
the Lucas- Washburn approach to include external pres- 
sure. As demonstrated, this allows a consistent analysis 
of resulting data. Such an analysis yields information on 
the Laplace pressure (if the pore radius is known, the 
contact angle then can be estimated) as well as the per- 
meability (and hence the slip length). The consistency 
of our description has been tested by simulations where 
all these quantities were obtained independently. Thus 
we have provided a framework which should be a useful 
guide for both experimental work on capillary filling and 
further related simulations. 
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